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In  the  last  six  months  we  were  able  to  complete  our  work  on  the  calculation 

of  the  macroscopic  optical  constants  of  an  assembly  of  clusters  of  random  orientation. 

As  mentioned  in  the  preceding  Reports,  a  cluster  of  spheres  is  a  model  scatterer  suitable 

to  approximate  the  properties  of  a  nonspherical  particle.  An  assembly  of  such  clusters 

should  therefore  be  a  good  approximation  of  an  actual  aerosol  composed  of  nonspherical 

particles.  The  work  done  in  the  last  semester  has  mainly  required  the  translation  of  the 

formulas  reported  in  the  paper  attached  to  the  preceding  Report  into  a  set  of  computer 
% 

programs.  As  this  has  been  done,  the  theory  has  been  applied  to  a  few  model  aerosols 

chosen  so  as  to  yield  a  general  in  sight  in  to  the  kind  of  the  results  obtainable  through 

our  approach.  A  careful  through  analysis  of  the  results  supports  the  reliability  both  of 

the  theory  and  of  the  computer  programs.  The  details  of  the  calculations  are  reported  in 

the  enclosed  paper^  which  has  been  submitted  to  "11  Nuovo  Cimento".  Our  results  have 

(2) 

also  been  presented  at  the  69. th  Annual  Meeting  of  the  Italian  Physical  Society  .  In 
the  subsequent  discussion  it  has  been  suggested  that  our  approach  could  be  applied  to 
the  study  of  the  absorption  from  clouds  of  interstellar  dust,  an  optical  medium  whose 
physical  properties  meet  the  requirements  for  the  full  applicability  of  the  theory. 

We  want  to  enfasize  that  the  programs  we  used  to  get  the  above  mentioned 
results  do  not  include  the  use  of  the  symmetry  properties  of  the  clusters.  Therefore  our 
next  step  shall  be  the  inclusion  of  the  symmetry  into  our  programs,  in  order  to  allow 
us  to  study  aerosols  composed  of  relatively  big  clusters.  The  above  mentioned  calcula¬ 
tions  are,  indeed,  the  most  complicated  that  can  be  done  without  the  help  of  group  theory. 
Furtermore,  we  notice  that  the  approach  we  used  to  connect  the  scattering  properties 
of  a  single  cluster  to  the  macroscopic  properties  of  the  aerosol  is  quite  elementary.  A 


3. 


possible  line  of  development  and  improvement  of  our  work  could  be  the  use  of  a  more, 
sofisticated  approach  for  the  connection  of  microscopic  and  macroscopic  properties. 
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Summary . 


^  A  method  to  calculate  the  macroscopic  optical  cons¬ 
tants  of  a  low  density  medium  consisting  of  a  cloud  of  identi¬ 
cal  nonsDherical  scatterers  is  presented.  The  scatterers  in 
the  medium  are  clusters  of  dielectric  spheres  and  the  electro¬ 
magnetic  field  scattered  by  each  of  the  clusters  is  obtained 
as  a  superposition  of  multipole  fields,  as  previously  proposed 
by  the  authors.  The  transformation  properties  of  the  spherical 
multipoles  under  rotation  allow  the  orientational-dependent 
terms  in  the  expression  for  the  forward-scattering  amplitude 
of  each  of  the  clusters  to  be  factored  out.  In  this  way  the  sum  of 
the  scattering  amplitudes  of  the  clusters  with  different  orien¬ 
tations,  needed  to  c  a  lc  ul  a  \  e  •  t  h  e  optical  response  of  the  medium,  ^ 


is.  greatly  facilitated  and  admits  a  'simple  analytic  expression 
in  the  case  of  randomly  oriented  clusters.  Results  of  calcula 


tions  of  the  optical  constants  for  a  few  model  media  are  prese 


nted  . 
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1  .  I NTRODUCTION 


It  is  well  known  that  the  key  quantity  for  the  evalu¬ 
ation  of  the  optical  properties  of  a  random  distribution  of 
scatterers  constituting  a  low  density  medium  ,  such  as  e.g. 
an  aerosol,  is  the  electromagnetic  scattering  amplitude  matrix 
of  the  individual  scatterers.  In  fact,  the  relevant  optical  pa- 
rameters  of  such  a  medium,  namely  the  refractive  index  and 
the  absorption  coefficient,  can  be  in  principle  obtained  from 


the  sum  of  the  forward-scattering  amplitudes  of  each  of  the 

(1  ,2) 


scatterers 


When  these  have  spherical  symmetry  their 

(1  ,2  , 


scattering  amplitudes  are  determined  from  the  Mie  theory 
30 


and  the  sum  leading  to  the  optical  parameters  presents  no 
additional  difficulty  on  account  of. the  assumed  low  densi¬ 
ty  of  the  distribution.  The  Mie  theory,  however,  does  not 
seem  satisfactorily  applicable  to  real  aerosols,  for  the  expe¬ 


rimental  data  support  the  conclusion  that  the  particles  of  the 

(4,5) 


most  common  aerosols  behave  as  nonspherical  scatteres 

In  this  paper  we  introduce  a  method  to  calculate  the 
relevant  optical  constants  of  a  medium  consisting  of  a  low  den 
sity  cloud  of  identical  nonspherical  scatterers  in  the  case 
that  each  of  them  is  a  cluster  of  dielectric  spheres.  These 


model  scatterers  proved  to  be  computationally  sound  and  to 
approximate  some  kind  of  nonspherical  scattering  centers: 


e.g.  they  hold  a  dependence  of  their  total  cross  section  on 
the  direction  of  the  propagation  vector  of  the'incident  field. 

(6  »  7 , 8 , 9)^  oreov  e  r  the  expression  of  the  scattering  amplitude 
of  each  cluster  can  be  given  in  terms  of  the  amplitudes  of  the 
multipole  fields  scattered  by  the  cluster  itself  and  does  not 
require  too  much  computational  effort.  The  literature  does 
not  report  calculations  of  the  optical  behaviour  of  a  medium 
of  nonspherical  scatterers,  but  in  the  case  they  are  all  ori¬ 
ented  alike^1  ^  .  I  ndecd  ,  identic  al  nonspherical  scatterers 
have  a  different  forward-scattering  amplitude  depending  on 
their  orientation.  This  greatly  complicates  the  sum  over  the 
scattering  amplitudes  needed  to  calculate  the  optical  response 
of  a  cloud  of  such  objects^\  The  medium  we  proposed  above 
allows  us  to  overcome  these  difficulties  and  to  calculate  its 
optical  constants  even  when  the  clusters  have  different  orien¬ 
tations  . 

Therefore.in  section  2,  we  recall  the  theory  to  cal¬ 
culate  the  amplitudes  of  the  multipole  fields  scattered  by  a 
single  cluster^’  ^and  put  the  final  equations  in  a  form  suita¬ 
ble  for  our  purposes.  In  section  3,  we  exploit  the  transforma 
tion  properties  of  the  multipole  fields  under  r  ot  a  t  i  on^ 1 1  ^a  nd 
succeed  in  factoring  out  of  the  expression  for  the  multipolar 
amplitudes  all  the  terms  depending  on  the  orientation  of  the 


cluster-  In  this  way  the  sum  of  the  scattering  amplitudes  of 
clusters  of  different  orientation  is  greatly  facilitated. 
Furthermore,  when  both  the  space  distribution  and  the  orienta¬ 
tion  of  the  (identical)  scatterers  are  random,  the  above  sum 
is  shown  to  assume  a  simple  analytic  expression.  At  last, in 
section4,  results  of  the  calculations  of  the  optical  constants 
for  a  few  model  media  are  presented.  In  particular  we  devised 
a  simple  model  medium  exhibiting  birefringence  and  dichroism 
even  when  its  clusters  are  randomly  oriented. 

We  conclude  these  introductory  remarks  observing 
that  no  dispersion  relation  is  taken  into  consideration  since 
our  treatment  is  mainly  concerned  with  the  effects  of  nonspheri 
city  of  the  scatterers.  For  the  sake- of  simplicity  neither  the 
case  of  the  clusters  of  conducting  spheres  is  included,  on  ac¬ 
count  that  such  an  extension  does  not  present  further  difficul¬ 
ties  or  advantages  for  our  purposes. 


2.  THE  FIELD  SCATTERED  BY  A  CLUSTER  OF 
SPHERES 

In  this  section  we  concisely  recall  how  the  mul  t  i  pol  a  r 
amplitudes  of  the  field  scattered  by  a  cluster  of  spheres  can 


5. 


(7) 

be  calculated,  referring  elsewhere  for  further  details 
In  the  following  both  the  spheres  and  the  matrix  including 
them  are  assumed  to  be  nonmagnetic,  dielectric,  homogeneo¬ 
us  and  i sotropic  . 

We  consider  a  circularly  polari  a  plane  wave  inci- 
ding  on  the  cluster,  i.e.  the  incident  fie  is 


li'^E  ( c  +-cne  )  e°-  -  , 
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with  n=±l  according  to  the  helicity.  For  our  purposes  we 
need  the  multipolar  expansion  of  the  field  of  eq.(2-l)  around 
the  center  of  the  a-th  sphere:  * 
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(2-2) 


In  the  above  equations  R^is  the  vector  coordinate  of  the  cen¬ 
ter  of  the  a-th  sphere  and  r  =r-R  ,  j,  are  the  Bessel  spherical 

r  -a  -  -  a  'l 

(2  11  12) 

functions  and  X  are  the  vector  spherical  harmonics  '  ’ 

-in 


6. 
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I  -/ 


We  expand  the  scattered  field  as 


*  .* 


'r\ 
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£S=  I  lA  „,„h,(kr  )x,  <f  )  +  B  Ivxh  (kr  )X  (f  )]  ,  (2-3a) 

aLM  ’  LM  L  0  ~LN  -a  n^LH"  L  a  "IK  *a 

<-?<nS=  I  [B  h  (kr  )X  (f  )+  A  ,  lv«h  (kr  )X  (f  )]  ,  (2-3b) 
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'while  within  the  a-th  sphere  the  field  ‘is 
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In  the  equations  above  A  ,  B  ,  ,  C  ,  „  and  D  ,  are  the 

n.oLM  n.oLM  n.oLM  n.aLM 

multipolar  amplitudes,  h^  are  the  Hankel  spherical  functions 

of  the  first  kind,  n  is  theCreal, constant)  refractive  index  of 

a 

the  o-th  sphere. 

(13) 

An  useful  addition  theorem  allows  us  to  expand 
the  field  of  eqs.(2-3)  in  the  vicinity  of  the  surface  of  the  a- 
th  sphere,  then  we  apply  the  tangential  continuity  conditions 
to  E  and  B  to  determine  the  multipolar  amplitudes.  .The  ampli¬ 
tudes  C  and  D  can  be  eliminated  to  obtain  the  system  of 
n,a  L  M  n,a  L  M 

linear  nonhomogeneous  equations 
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In  eq s  . (2-4) 
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where  b  is  the  radius  of  the  a-th  sphere  and  j  =j  (kr  );h  = 
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F  urther more 
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where 
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and  the  C' s  are  the  Clebsch-Gordan  coefficients 
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are  the  Gaunt  integrals 

The  solutions  of  the  system  (2-4)  are 
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In  the  equations  above 


r(A)  p  ik.R  ,  ,-(**)  T  (  ab  ) 

=  -  )  e  ~  -a'  ( T  +  n  T  ) 

n,alH;l'N'  L  alHjo'l'H'  al*;a'l'K' 


-6b) 


-7) 


I  n  e 


8a) 


8b) 


» 


9. 


j(  e) 

n,oLN;L'M 


(  =  - 

.  ik-R 
I  e  -  - 

a* 

the 

matrix 

( AA) 

t(AB) 

(BA) 

•t(bb) 

(B A )  +  (BB)  J 

aLH;  a'L'N*  n  aLH;o'L'H' 


is  the  inverse  to  the  matrix  of  the  system  (2-4) 
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obviously  all  the  submatrices  have  the  same  number  of  rows. 
For  our  purposes  we  need  the  multipolar  expansion  of  the 
scattered  field  around  the  origin  of  the  axes.  At  a  distance 
larger  than  the  overall  size  of  the  cluster  the  addition  theo¬ 
rem  already  used  above  allows  us  to  write 
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In  eqs.(2-9)  then’s  and  the  X'  s  are  given  by  eq.(2-6a)  and 
(2-6b),  respectively,  with  substituted  for  and  the  omis¬ 
sion  of  the  factor  (  1-i  ,)  within  the  G's  of  eq.(2-7).  With  the 

oa'  ^ 

help  of  eqs.(2-8),  eqs.(2-9)  become 
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Concluding  this  section  we  specifically  enlighten 


two  points.  The  first  remark  regards  the  variables  on  which 
the  U's  and  W's  in  eqs.(2-10)  depend.  The  definition  (2-2) 
show  that  the  W's  depend  only  on  e  and  $  ,  the  polar  angles 


of  k.  Tracing  through  eqs.(2-ll)  and  related  definitions, 
one  finds  that  the  U's  depend  only  on  the  magnitude  of  k  ,  the 
n^'s,  the  b^'s,  and  on  the  components  of  the  R^'s  .  Then  by 
eqs.(2-10)  we  succeeded  in  expanding  the  multi  polar  amplitu¬ 
des  as  linear  combinations  of  the  geometrical  factors  U  with 
the  orientation  -  dependent  coefficients  W.  The  second  remark 
regards  the  order  of  the  system(2-4).  We  have  2NL..CL..+  2) 
equations,  where  N  is  the  number  of  the  spheres  in  the  cluster 
and  L  is  the  maximum  L  to  be  retained  in  the  multipolar  expa- 
ns  ions.  Thus  the  order  becomes  easily  very  large  and  some 
advices  should  be  applied.  For  example  one  could  exploit 
the  symmetry  of  the  cluster,  as  we  reported  elsewhere^ 


3.  THE  REFRACTIVE  INDEX  AND  THE  ABSORPTION 
COEFFICIENT  OF  THE  MEDIUM 
To  begin  this  section  we  make  a  remark  about  the 
field  inciding  on  each  cluster  of  the  medium  we  are  conside¬ 
ring:  because  of  the  multiple  scattering  processes  this  field 
should  also  include  the  fields  scattered  by  all  the  other  clu¬ 
sters.  However,  disregarding  them  amounts  to  a  fairly  acce¬ 


ptable  approximation  when  the  concentration  of  the  scatterers 


is  sufficiently  low  .  Since  we  deal  just  with  this  low  den¬ 

sity  limit,  eqs.(2-10)  give  also  the  multipolar  amplitudes  of 
each  cluster  belonging  to  the  medium.  Nevertheless ,  nobody 
could  actually  calculate  the  scattering  amplitudes  of  all  the 
clusters,  one  at  time,  referred  to  an  unique  system  of  coor¬ 
dinates,  and  sum  them  to  get  the  optical  response  of  the  me- 
'dium.  This  formidable  difficulty  can  b'e  overcome  in  two  steps 
however,  through  the  use  of  the  assumed  identity  of  the  clu¬ 
sters  and  of  the  transformation  properties  of  the  multipole 
fields  under  rotation. 

Fist  of  all  we  define  a  system  of  coordinates  l  to 
which  the  whole  medium  can  be  referred.  Furthermore,  let 

-fv) 

us  refer  each  cluster  to  a  local  set  of  axes,  l  for  the  v-th 
cluster,  chosen  so  that  the  local  coordinates  of  the  centers 
of  the  corresponding  spheres  in  different  clusters  are  identi¬ 
cal.  We  also  associate  to  the  system  of  coordinates  J  ^ 

with  the  same  origin,  but  with  the  axes  oriented  as  the  axes 
of  l  .  Now  we  can  write  the  electric  field  incident  on  and 

scattered  by  the  v-th  cluster  with  reference  to  the  system 

f  (v) 

1  as 


E  W(vMj  (k?(v))X  MCr(v))  + 

O  n  t  k  n  l  n 


+  n— 7xj  (kf(v))  X  (r(v))] 


(3-la) 


13. 


=  1  [jl h/k?(V^x. /r(v))  +  C(v)  ^7xh  (k?{v))x  (?(vV]  ,  (3-1  b) 
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where,  according  to  eqs(2-10),  the  amplitudes  are 


/  M  =  £  ik.Rvj  (A)  (v) 

C  ■sin  ^  L  -  I  M  .  II  Ml  -  II 


n.  L'H* 


(3-2a) 


C  lv)  =  E  e^-?vTu(B)  W(v) 
n.L«  O  *  n, LM  ;L'M'  n.L'H* 


(3-2  b) 


-  (v) 

Since  the  components  of  R#  are  independent  of  v,  neither  the 

U‘s  in  eqs.(3-2)  depend  on  it;  on  the  contrary  the  W's  depend 

_  f\)  ~(v)  i-k  •  R 

on  v  through  8fc  and  .  The  phase  factor  e  -  ~v  refers 

mm 

the  amplitude  Eq  of  the  plane  wave  in  the  origin  of  J  to  the 

"(v) 

amplitude  in  the  origin  of  J  (  of  vector  coordinate  R^  )  .  The 

electric  fields  can  also  be  written  directly  with  reference  to 
t  <v) 

I  as 


E(°  =  E  etv’BvJw  [j  (kr(v,)X  (r(v))  + 
~ n  v  o  ‘  l*  l  -m  - 


+  n-^v  x  j  (kr(v))X  (r(v,)J 

K  L  ~  l  A  ~ 


C3-3a) 


It  is  to  be  noticed  that  the  W’s  in  eq.(3-3a)  do  not  depend  on 

v,  viz.  on  the  orientation  of  the  cluster;  they  actually  depend 

(v) 

only  on  the  polar  angles  of  k  with  respect  to  the  axes  of  l 
which  by  definition  are  parallel  to  the  axes  of  l  .  Moreover, 
on  account  of  their  transformation  properties  under  rotation, 

the  vector  spherical  harmonics  in  eqs.(3-l)  and  (3-3)  can  be 

,  J  .  .  .  (11 ,12) 

referred  to  each  other  by 


X  (?{v))  =  l  DU)(e)  X  (r  {v)) 

-I  «  ~  L  M*M  v  — L  M*  — 


(3-4) 


where  the  argument  e  of  the  representation  matrices  of  the 

full  rotational  group  represents  the  Euler  angles  of  the 

”(v)  (v) 

rotation  transforming  l  into  l  .  Thus  ,  by  substituting  eq 
(3-4)  into  eqs.(3-l)  and  comparing  with  eqs.(3-3)t  we  get 


w  .a,*  I  DU)  (e  )  w*v) 

n,Uf'  *  v  n.L« 


whose  inverse  i  s 


W  (  *  «  l  W  D(^(e  ) 
n.ut  4  n.LN'  v 


(3-5) 


i  d5!c  >  ^v! 


(3-6a) 


(3-6b) 


£(v)  =  l  o(Lle  )  ?v) 

n.w  J  m<n  VV  n.iN 


Ultimately,  eq s  .  (3-6) ,  (3-2)  and  (3-5)  yield 


9{  ]  -  E  e  -  l  l  U(  1  — W  -  D  L  (e  )  D(L)-(e  ) 
H.LB  O  --  L  n,LN';L*  n.LN"  K  H1  v  H"H  v 

U  »'«" 


(3-7  a) 


and 


B  (v)  *  E  e‘C-*-vI  l  U  (B)  — W  -  D{L)(e)  D(^*(e  ) 

n,L«  o  i-  ‘Nn  n.LttU"  n.tN"  NM*  v  >fV 


(3-7  b) 


Eqs.(3-7)  evidence  that  the  multipolar  amplitudes  of  the  scat¬ 
tered  field,  when  referred  to  I  ^  ,  depend  on  the  orientation 
of  the  cluster  only  through  the  elements  of  the  matrices  tfL\e^). 
The. other  quantities  in  eqs.(3-7),  with  the  exception  of  the 
phase  factor,  do  not  depend  on  v  and  can  thus  be  calculated 
once  for  all.  This  result  yields  great  advantage  for  the  calcu¬ 
lation  of  the  optical  constants  of  the  medium.  Indeed,  the  re¬ 
sults  obtained  thus  far,  eqs.(3-7),  can  be  related  to  macro¬ 
scopic  optics  in  a  way  that,  although  elementary,  will  maintain 
fairly  visible  the  effects  induced  by  the  lack  of  spherical 
symmetry  of  the  scattering  centers. 

The  normalized  f  or  w  a  r  d  -  s  c  a  1 1  e  r  i  n  g  amplitude  of  the 


v-th  cluster  is 


(2) 


(3-8) 


0  - - —I  (-<.)  X  (R)  +L  RxX  (R)]  . 

v,n  evS*  at*  n*L"  “LH  “  n*L"  “  ~L"  - 
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Defining  f  =u**f  ,  where  u*  =  -  is  the  (complex)  unit 

v,n  n  ~n'  ~v,  n  -  n'  j-  v 

polarization  vector,  eqs.(3-S)  and  (2-2)  yield 


v.n'n 


-  e'^.^v  Jw*  (^{v)  +n'/?(v)  ) 

8*  -ckE  L  N  T'',Ln  n,LH  T),LM 


(3-9) 


These  quantities  are  in  turn  related  to  the(complex)  refracti 
ve  index  of  the  medium  by 


?  -  =  «  ,  +  —  If  , 

n'n  n'n  v.n'n 

rV  v 


(3-10) 


where  the  sum  includes  the  clusters  within  the  volume  V.  We 

recall  that  n  = /?e.(  7?  )  and  y  =  2k Jm(7l  )  are  -the  macroscopic 

n  nn  n  nn  r 

refractive  index  and  the  absorption  coefficient  of  the  medium 

respectively;  the  off-diagonal  elements,  77  ,  account  for  the 

nn' 

rotation  of  the  plane  of  polarization.  Since  f  =  f  .  when 

r  v.n'n  v.nn 

the  identical  clusters  v  and  v'  have  also  the  same  orientation 
eq.(3-10)  can  be  rewritten  as 


I  N.  f 


(3-11) 


where  N.  is  the  number  density  cf.  clusters  having  the  i-th 
orientation.  Provided  the  low  density  of  the  medium  is  not  too 
low,  eq.(3-ll)  can  be  replaced  by 


6  +  N(e)  f  ,  (e)  d  e 

nn  . .  nn 

k* 


(3-12) 


This  results  in  a  further  substantial  simplification  in  the  case 

that  also  the  orientation  of  the  clusters  is  random.'  In  this 
N 

case  N(e)  =  -  ,  where  N  is  the  number  density  of  all  the  clu- 

8  IT  * 

sters ,  and  we  have 


d  (L  it  H';  l  M'M)  =  N 


2U1 


6-6-6-  , 
L  L  MM'  M'M 


(3-15) 


on  account  of  the  orthogonality  properties  of  the  rotation  ma- 

(0(11) 


trices  D  .  Then  eqs.(3-13)  with  eqs.(3-7)  and  (3-15)  be¬ 


come 


N  I  u(A) . W 


H.LM  2L  ♦!  n.LMUM'  n ,  L  M 


B 


N  iu"»  w 


n,  LM  2L  ♦!  n . LM •; L M 1  n.LM 


The  limiting  case  when  all  the  clusters  have  the  same  orienta¬ 
tion  is  trivial.  Assuming  whe  have  N(e)  =  N  fi(e)  and 

eq  .(3-1 4)  give 


d(l  M  MU  M'M  )  =  ND-L)-(0)D(L)(0)  =  N  6-  -  6 

M'M  MM'  M'M  MU' 


(3-16) 


Using  eqs.(3-13),  (3-7), (3-16)  and  (3-2),  we  'ind 


where  e  -  ~v  cancels  out  with  the  phase  factor  contained  in 


VV)  and  Z?(v) 
n.LH  n.LH 


4.  RESULTS 

We  applied  the  theory  developed  in  the  preceding  sec¬ 
tions  to  a  few  model  media  composed  of  clusters  of  rather  sim¬ 
ple  structure  with  random  space  and  orientational  distribution. 
Each  medium  has  been  considered  also  for  like  orientation  of 
clusters  with  several  directions  of  k,  the  wavevector  of'  the  in¬ 
cident  field.  This  allows  us  to  show  the  effect  of  the  orienta¬ 
tional  disorder  of  the  scatterers  on  the  optical  constants  of 
the  medium.  Table  I  reports  the  parameters  defining  the  clu¬ 
sters  of  each  medium  individuated  by  a  conventional  name  to 
which  we  will  refer  afterwards.  Of  course  the  names  are  descri 
ptive;  so  "Dichro"  denotes  a  simple  medium  capable  to  show 
an  evident  birefringence  and  dichroism.  The  coordinates  of 
the  centers  and  the  radii  of  the  spheres  as  well  as  1 /k  are 


measured  in  arbitrary  units.  The  radius  of  one  of  the  spheres 


has  been  labelled  as  b  and  chosen  to  enter  the  parameter  x  =  kb 
versus  which  all  the  results  will  be  reported.  Though  x  could 
clearly  be  related  to  the  ratio  between  the  size  cf  the  cluster 
and  the  wavelenght  of  the  incident  f  i  e  1  d ,  f  u  r  t  h  e  r  information 
on  this  point  can  be  gained  by  the  parameter  C=kbo»  where  b^, 
the  radius. of  the  smallest  sphere  containing  just  one  cluster, 
is  also  reported  in  Table  1. 

The  results  of  our  calculations  are  better  evidenced 
introducing  the  quantities 

n  -  1 
n 

P  n  (°) 

n  -  1 

where 


„(o)  (o)  .  Co). 

7}  =  (  n  +  cif  )  6  , 

n  n  nn 


are  the  elements  of  the  refractive  index  of  the  medium  that  one 
obtains  by  the  Mie  theory  when  the  spheres  in  each  cluster 
are  considered  as  independent  scatterers,  i.e.  disregarding 

themultiple  scattexing  processes  even  among  the  spheres  inthesame 

(7) 

cluster  .  p  and  g  are,  in  a  sense,  normalized  quantities 
n  n 


and  do  not  depend  on  any  concentration  of  scatterers  we  choose 
in  the  low  density  range.  All  our  calculations  were  performed 
for  x  ranging  from  0.001  to  1.0  but  all  the  results  are  actually 
reported  for  x  between  0.1  and  1.0,  the  most  interesting  ran¬ 
ge  as  far  as  the  effects  of  orientational  order  and  disorder 
are  concerned.  It  is  worth  noticing,  in  this  connection,  that 
when  x=1.0,  C  reaches  2.  for  ’'Oxygen'*,  3*5  for  "Water"  and 
3.1  for  "Dichro".  Therefore  the  wavelenght  is  always  larger 
than  the  radius  of  the  spheres  but  becomes  comparable  to  the 
overall  size  of  the  clusters. 

Figs.  1  and  2  report  p^  and  g^,  respectively,,  for 
"Oxygen".  The  solid  curves  refer  to  a  medium  whose  clusters 
are  oiriented  alike  for  several  values  of  the  polar  angle  , 
which,  in  this  case,  gives  the  angle  between  k  and  the  axes 
through  the  centers  of  the  spheres.  Quite  identical  curves 
could  be  drawn  for  p  ^  and  g  As  we  expected  the  dependen¬ 
ce  of  p  and  g  on  the  orientation  of  the  clusters  increases 
n  n 

with  x(5).  The  broken  curves  in  figs.l  and  2  give  pi  and  g  , 
respectively  ,  when  the  clusters  are  randomly  oriented;  they 
carry  no  label  because  for  random  orientation  the  results 
are  independent  of  the  direction  of  k.  It  is  worth  noticing  the 
crossing  of  all  the  curves  in  fig.l,  including  the  broken  one. 


at  a  single  value  of  x.  This  result  suggests  that,  as  far  asp 

■  n 

is  concerned,  there  exist  a  wavelenght  at  whtch  the  medium 

behaves  as  a  cloud  of  spherical  scatterers.  This  feature  is 

shared  in  common  with  other  media  and  may  occur  also  for  the 

curves  of  g  . 

n 

Figs.  3  and  4  report  p^  and  g^,  respectively,  for 
"Water".  Solid  and  broken  lines  have  the  same  meaning  as  befo¬ 
re  and  their  behaviour  suggests  analogous  considerations.  The 

solid  curves  are  labelled  by  the  value  of  6,  ,  while  ♦  =90°  . 

k  k 

Of  course  the  calculations  were  actually  performed  for  seve¬ 
ral  values  of  We  also  s  t  re  s  s  that  the  crossing  of  all  the  cur- 

ves  occurs  both  for  p1  and  g^  though  at  different  values  of  *. 

The  last  medium  we  consider  in  this  paper  is  "Dichro 
whose  Pj  and  gj  are  reported  in  figs. 5  and  6,  respectiv  ely . 

The  complete  lack  of  rotational  symmetry  of  the  clusters  in 
this  medium  prevents  all  the  curves  to  cross  at  a  single  point 
but  all  the  crossings  occur  within  a  small  range  of  x.  The  most 
striking  features  of  this  medium  are  its  birefringence  and  di¬ 
chro  is  m  which,  though  strongly  attenuated,  do  not  vanish  even 
when  the  clusters  are  randomly  oriented.  We  report  in  figs.  7 
and  8  the  quantities 


respectively,  for  n  =  l;  clearly  6  r-B  and  d  =  -d  .  The  differ- 

—  n  n  —  n  n 

rent  scale  used  in  figs. 7  and  8  for  the  case  of  random  orienta¬ 
tion  should  be  noted. 

The  off-diagonal  elements  of  7}  ,  eq.(3-10),  were  also 
calculated  for  all  the  cases  reported  above  and'found  to  be  non¬ 
vanishing  but  when  expected  according  to  the  discussion  of 
Newton^^.  These  results  together  with  those  reported  above, 
support  the  reliability  of  the  theory  developed  in  this  paper 
to  approximate  the  optical  behaviour  of  a  . cloud  of  nonspherical 
model  scatterers.  We  stress  that  our  theory  alsoapplies  when  the 
orientational  distribution  of  the  clusters  is  not  completely  ra¬ 
ndom.  In  this  case  it  is  necessary  to  estimate  the  d-integrals, 

eq.(3-14).  In  particular,  if  N(e)  is  proportional  to  a  generali- 

(14) 

zed  spherical  harmonic  the  d-integrals  can  still  be  given 

analitic  ally .  In  any  case,  nothing  prevents  their  use  as  dispo¬ 
sable  parameters,  thus  transforming  the  theory  into  a  semiem- 


pirical  method  . 
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R  i  a  s  s  u  n  t  o 


Viene  presentato  un  metodo  per  il  calcolo  delle  costan 

ti  ottiche  mac  roscopic  he  di  un  mezzo  di  bassa  densita  costitui- 

to  da  una  nube  di  diffusori  identici  non  sferici.  I  diffusori 

del  mezzo  sono  "cluster"  di  sfere  dielettriche  ed  il  campo  elet- 

* 

tromagnetico  diffuso  da  ciascuno  di  essi  e  ottenuto  come  sovrap- 
posizione  di  campi  multipolari  ,  come  precedentemente  proposto 
dagli  autori .  Facendo  uso  delle  proprieta  di  trasformazione 
dei  multipoli  sferici  sotto  rotazione,  i  termini  dipendenti  dalla 
orientazione  nell'espressione  d  e  1 1 '  a  m  pi  e  z  z  a  di  diffusione  in 
avanti  vengono  f  att  ori  z  z  ati  .  In  questo  modo  la  somma  delle  am- 
piezze  di  di'f  f  u  si  on  e  di"cluster"  con  differenti  ori  ent  a  zi  oni  , 
necessaria  per  il  calcolo  della  risposta  ottica  del  mezzo,  vie¬ 
ne  molto  facilitata  e  assume  una  semplice  espressione  analiti- 
ca  quando  i  "cluster"  sono  orientati  a  caso.  Vengono  anche 
presentati  i  risultati  del  calcolo  delle  costanti  ottiche  per 
qualche  mezzo  modello. 
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